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Abstract 

While the investors' responses to price changes and their price forecasts are well accepted 
major factors contributing to large price fluctuations in financial markets, our study shows 
that investors' heterogeneous and dynamic risk aversion (DRA) preferences may play a more 
critical role in the dynamics of asset price fluctuations. We propose and study a model of an 
artificial stock market consisting of heterogeneous agents with DRA, and we find that DRA 
is the main driving force for excess price fluctuations and the associated volatility clustering. 
We employ a popular power utility function, U(c, 7) = c with agent specific and time- 

dependent risk aversion index, 7i(t), and we derive an approximate formula for the demand 
function and aggregate price setting equation. The dynamics of each agent's risk aversion 
index, ji(t) (i=l,2,...,N), is modeled by a bounded random walk with a constant variance 
5 2 . We show numerically that our model reproduces most of the "stylized" facts observed in 
the real data, suggesting that dynamic risk aversion is a key mechanism for the emergence of 
these stylized facts. 

1 Introduction 

There have been many attempts to construct models of financial markets and to understand the 
key statistical features of financial time series. It remains a great challenge, however, due the 
inherent complexity of the financial market, to develop a parsimonious market model that can 
reproduce all key "stylized" facts observed in real financial data and provide insights into the 
market mechanism for the emergence of these stylized facts. One of the promising approaches 
is agent-based modeling, which has reproduced and explained the emergence of some of the 
"stylized" facts. The agent-based modeling provides an ideal framework for investigating the 
impact of investors' behaviors on the price dynamics from many different perspectives; it has 
become an indispensable tool for understanding the price dynamics of financial markets 1 . With 
agent-based modeling, one can model and investigate, for example, how investors make their price 
forecasts and how their price forecasts influence the price fluctuation (Arthur, Holland, LeBaron, 
Palmer, and Tayler (1997), DeLong, Shleifer, Summers and Waldmann(1990), Levy, Levy and 
Solomon (2000)), how investors respond to price change (Lux and Marchesi(1999), Caldarelli, 
Marsli and Zhang (1997)), and how investors form and change their market beliefs (Brock and 
LeBaron (1996), Barberis, Shleifer and Vishny (1998)). In this paper we use agent-based models 
to study how investors' fluctuating risk preferences affect the price dynamics. This important 
issue has not been fully explored. 

Financial markets present many important and challenging problems. First, the market 
consists of intelligent, competing and heterogeneous agents who, with different beliefs in the 
market, different abilities to acquire and process market information, and mutually conflicting 
interests, try to make investment decision for their own benefits. Second, each agent's decision 
depends on his estimate of price expectations of other agents who also make their own estimates; 
this precludes expectations being formed by deductive means and leaves inductive reasoning as 



1 LeBaron(2005) is a good recent survey on the field 
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the only choice (Arthur et al 1997). A market of agents employing inductive reasoning often 
exhibits irrational herding behavior (Bak, Paczuski and Shubik (1997), Cont and Bouchaud 
(2001)), resulting in excessive price fluctuations or sometimes market bubbles and crashes. In 
such market agents' sentiment and degrees of risk aversion play a critical role in determining 
its price dynamics. Third, agents can learn and therefore adapt their strategies dynamically to 
improve their performance; this exacerbates the unpredictability of the markets. The change 
in an investor's strategy or behavior can be the cause or the result of the investor's changing 
sentiment (Barberis and Shleifer 1998), represented by (pessimistic) under-reaction or (optimistic) 
over-reaction to market dynamics driven by arrival of new information and changing market 
macro/micro-environment. These characteristics of financial market (competing with different 
beliefs and conflicting interests, interdependence of price expectations, and unpredictable changes 
of risk aversion) may lead to the formation of so-called noise traders behaviors. DeLong et al 
(1990) have found this a (behavioral) source for price to diverge significantly from the fundamental 
value in financial markets - the so called "noise trading effect". A good model, needless to 
say, must successfully address these challenges; more importantly, it must be able to produce 
simulation time series that can capture the key "stylized" empirical facts observed in real financial 
time series. This paper reports our efforts in constructing such a model. The key ingredient of 
our model is the inclusion of investors' changing degrees of risk aversion. We show that this is 
the main cause of excess stock price fluctuations and the associated volatility clustering. 

We outline here our model of interacting heterogeneous agents. We first consider a baseline 
model, in which the agents use past price information to form their sets of future price expec- 
tations. The agents are adaptive as their price expectations are not based on one particular 
estimator but are determined by their best-performance estimators which may change from time 
to time. The agents also use their erroneous stochastic beliefs (DeLong et al. (1990)) in the market 
to make price adjustment on their best forecasts, which are assumed to be normally distributed. 

In our baseline model, we assume that the agents have a decreasing absolute, but constant 
relative, risk aversion (DARA and CRRA) utility function, U(c) = (c 1_7 — 1)/(1 — 7). The 
existing agent based models of stock market, such as the Santa Fe artificial stock market model, 
typically use a constant absolute risk aversion (CARA) utility function, U{c) = — e _Ac , for the 
price setting equation can be easily derived under such utility. As the focus of our paper is on risk 
aversion, we choose to use the well accepted (DARA) power utility function. Although a simple 
analytic formula for the demand function under this utility is not available, we have derived a 
general functional form and a rather good approximation for the demand function. Like the 
SFI market model, our baseline model market with a choice of the parameters corresponding to 
normal market conditions, exhibits some excess volatility, but not to the extent of the volatility 
observed in real markets. In addition, there is little enhancement of volatility clustering at high 
volatility regime, which is observed in real market data (Chen, Jayaprakash and Yuan (2005)). 

By simply allowing investors to change their risk aversion attitudes, we obtain excess volatility 
and volatility clustering in very good agreement with real market data. The implication is clear: 
dynamic risk aversion (instead of fixed constant risk aversion) is directly responsible for excess 
volatility and the associated clustering. Specifically, in our DRA model, which is built from our 
baseline model, all agents have the power utility functions ((c 1-7 — 1)/(1 — 7)) but with different 
and time varying risk aversion indices (degrees), 7^, which we assume to follow an independent 
bounded random walk with a variance e> 2 . We will show that the magnitude of excess volatility 
is directly related to 5 2 . With such DRA our model market exhibits most of the important 
statistical characterization of real financial data, such as the "stylized" facts related to excess 
volatility (Mandelbrot (1963), Fama (1963), Bouchaud and Potters(2000), Mantegna and Stanley 
(1999), Cont (2001)) and volatility clustering (Mandelbrot (1963), Fama (1965), Engle (1982), 
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Baillie, Bollerslev, and Mikkelsen (1996), Chou (1988), Schwert (1989), Poterba and Summers 
(1986), Chen et al. (2005)). 

We have also studied the impacts of the dynamic risk aversion on the market dynamics using 
a few other baseline models, including the SFI market model, and we found similar results. This 
suggests that our results on dynamic risk aversion are rather generic. 

The next section contains a derivation of the price equation under the power utility function 
for risk aversion. Section III describes our baseline model with a fixed constant risk aversion. 
Section IV introduces our DRA model. Section V reports the results from numerical simulations 
of our model. Section VI considers a DRA model built with the SFI market model as the baseline 
model. The last section summarizes. 



2 Demand Function and Price Setting Under the Power Utility 
Function 

We consider a market of N heterogeneous agents who form their subjective expectations induc- 
tively and independently based on their investment strategies. There are two assets, a risky stock 
paying a stochastic dividend with a limited supply of N shares 2 , and a risk-free bond paying a 
constant interest rate, r, with infinite supply. All agents have the same form of power utility 

function, U(ct;ji t t) = *i_ T - , but they have their own time-dependent risk aversion indices 
denoted by 7^. At each time step t, every agent decides how to allocate his wealth between the 
risk-free bond and the risky stock. Since the values for both the dividend payment and the stock 
price at the next period t + 1 are unknown random variables, the investors can only estimate the 
probability of various outcomes. Assume each agent's estimation at time t of the next step's price 
and dividend is normally distributed with the (conditional) mean and variance, En\p t+ i + dt+i] 
and crf t {i = 1,2, ...iV) respectively. It can be shown, by optimizing the total utility, that the 
demand of agent i for holding the share of the risky stock is approximately, 

where pt is the stock price at time t, 7^ is agent i's index (degree) of risk aversion, and af t the 
conditional variance of price estimation. 

The market clearing condition: Di y t+ T = Tld Dij = N can be used to determine the 
current market price and relate the price at time t + T, pt+ T to the price at time t, pt' 

E iit \pt+i + d t +i] 

Pt = ,r , ^ 1 • ^ 



and 



spN Ei, t+T \p t+T+1 +dt +T+1 ] _ y^TV 1 

^ 7i,t+T<r?, t+T (l+r) ^ 7i,t+Ttr?, t+T . . 

Pt+T y^jy gj.tbt+i+dt+i] y^iV 1 Pt ' (0) 



It can be seen from the demand and price equations that the degree of the agent's risk aversion 
plays an important role. 



2 In practice the number of shares is never the same as the number of agents; here we set the two numbers the 
same for the sake of convenience and setting them different does not change the results 
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We now show the derivation of the above equations. Assume at time t agent i's consumption 
is Cj t, and he invests a portion x of his current consumption in the risky asset. His total utility 
function defined over the current and future values of consumption is: 

U(c i t,c i t+i,'Yi,t) = U(ci t ,y it ) +U( Ci ' tJ £ X ,j it ), (4) 

where the consumption at time t + 1 can be written as 

c i)t +i{x) = Cj it [(l - x)R f + xRi tt+ i\. (5) 

Here Rf = 1 + r is the gross risk- free return and R% t+i the gross return on the risky asset. Agent 
i determines the amount of his investment on the risky asset, x, by maximizing his total utility, 
Eqn (|4"|). The maximization problem can be written as 

max Etpia^t) + C/(^M, 7M )] = max^p(^, %t )}, (6) 

The last equality follows because the utility U (c^t, 7i t) is known at time t and it does not contain 
x. 

The power utility function is given by, 

cjr M - 1 

tf(ci, t ;7i,t =-r ( 7 ) 

1 - 7t,t 

Substituting Eqn Q into Eqn (JJJ) and then inserting it back to Eqn Q, the maximization is now 
given by 

c 1-7M 

max _W ^{[1 + f <itf ix] 1 "^.*} (8) 

21 1 - 7i,i 



where fj j+i = *' "U £ is the present (discounted) value of the net return at next time step 

t+1. Assuming that agent i's prediction errors of fa+i are (conditionally) normally distributed: 



f i)t+ i = Et(f i)t +i) + Zi )t = r e it + z ijt (9) 

where rf i = Etftn+i) is the conditional expected net return, at time t, of the next time step; 
and the error of estimation is z% t ~ N(0, crf t ). The maximization becomes: 

max^{(l + xr| f + xz it ) 1 ~ 7 '- t } = max{ / e 2<Ti .* 2, -y i:t )dzij}, (10) 

where f(zi t t; x, 7^) = (1 + xr? t + sc8i,t) 7<,t . In writing down the above equation we implicitly 
assumed (1 + 2T? t + 22j t) > 0, which is a necessary requirement for the power function to be a 
valid utility measure. The above integral can be further simplified to 

1 f°° 2 

— / e~ Zi >* f{y/2ai,tZi,f,x,"1i,t)dzi 7 t (11) 
V 71 " i-00 

This Gaussian integral cannot be evaluated analytically, but it can be approximated by an ex- 
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pansion based on the roots of Hermite Polynomial H^ n \^) as: 

/oo n 
e- z ^f{y/2a i>t Zi,t;x, li>t )dz i)t = ^ X^> f (V2a ijt ^ ; x, %t ) (12) 
-CO 7. -i 



k=l 



where = 1,2, ...,n) are the coefficients of the summation, and £^ are the roots of nth 

Hermite polynomial H^ n \^). Performing the maximization by setting the derivative with respect 
to x equal to zero, we obtain the following equation: 



X>J°[1 + (rf, + V2a i;t ^)x}-^ x (rf >t + v^tfj ) = ( 13 ) 

k=l 



Since |rf t | <C 1, £^ ~ ^(1) 3 an d <^i,t ^ 1 f° r the typical time step of one day or shorter, the 
above can be approximated as 



£ A< n) [1 - 7M (rk + v^<,4 n> )*] >< K* + ^r*,4 B) l = 0- ( 14 ) 



fc=l 

(2) (2) / — (2) 

Here we consider an approximation with n = 2. Note that A^ = A 2 = A, y2q = 
— \/2^2 ( = -Qj the optimal demand of agent i of risky stock can then be obtained as 

n = r t,t „ r i,t E i)t \pt+i + dt+i] - (l+r) Pt 

i ' t ^ 7mI(^ 2 + <J~7m^ 7M^*(l + r)ft ' 1 j 

which is the Eqn. ([T]l. In writing down the above approximation we assume (rf t ) 2 <C cr? t , which 
is certainly true when the time step is one day or shorter. 

To get a more accurate approximation of the demand function, higher order Hermite polyno- 
mial roots are needed in the summation approximation to the integral in Eqn. (|11|), It can be 
shown that, the demand function Xi in a higher order approximation is exactly the same as Eqn. 
(|15j). except for an overall constant factor. 



3 The Baseline Model 
3.1 Price prediction 

To use the demand and price setting function derived in the previous section, one still need to 
incorporate each agent's prediction of the payoff at the next time step, Et(pt+i + dt+i)- We 
assume all agents use the past price information for price forecasting. The simplest way is to 
calculate a moving average of the available past prices and use it as a proxy of price forecasting 
for the next time step t+1. Since investors may have different investment horizons and evaluation 
strategies, they may use different time lags for their calculation of the moving average of past 
prices, implying that they have heterogeneous memory lengths (Levy et al. (1994)). We consider 
each agent has his own M sets of predictors so that he can choose the best one for forecasting the 
price at the next time step. Each price predictor Eij(p t +i + dt+i), (i = 1, 2, N; j = 1, 2, M) 
is made of a moving average of past Lij prices with a subjective erroneous stochastic adjustment: 

E ld ,t(pt+i + dm) = MA^ t = MAi iiiW (l - -!-) + — (pt + d t ) + eij (16) 
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where e^j ~ iV(0, <J p +d) is Gaussian random variable. 

The conditional variance of the estimation, af t , is assumed to update with a moving average 
of the squared forecast error: 

a\ ht = (1 - 0)^-1 + 0[(pt + dt) - £ lli)t _i{pt + d t )} 2 , (17) 

where #(0<#<Sl)isa weighting constant. 



3.2 Dividend process 

The dividend process is assumed to be a random walk: 

d t = d t -i +r d + e t , (18) 

where et is an i.i.d. Gaussian with zero mean and variance cr^; is the average dividend growth 
rate. Note that the dividend process in a real stock market may be more complicated than what 
we assumed here and it may vary from stock to stock. But our results are not sensitive to the 
choice of a dividend process. 



3.3 The price setting equation of the baseline model 

For heterogeneous agents with fixed constant risk aversion, the demand function Eqn. Q and 
price setting Eqn. @ and © can be written as: 

A t = ^ [pm+ / m] " Pf(r + 1) (19) 
Wi,t\ 1 + r )Pt 

and 

E i>t \p t+1 + dt+i] 



Pt = N + T N ! (2°) 



v^iV 1 / Ei, t+T [p t+T+ i+d t+T+ i] ^ 1 \ 

It can be see that the risk aversion indices of the agents play the role of weighting factors in 
the price setting equations. 

If the agents are homogeneous in risk aversion, the above can be further simplified to: 

E ht [p t+1 + d t+ i] -p t (r + 1) 

wfA 1 + r )Pt 

and 



From the above equations we see that the risk aversion index 7 only affects the overall level 
of demand but not its fluctuations. It also does not contribute to the price fluctuations (between 
the time t and time t + r). Thus in the case of homogeneous and constant risk averse agents, the 
main source of the price fluctuations is from the investors' price forecasting. The baseline model 
does not incorporate investor's changing sentiment. As a consequence, the price fluctuation is 
expected to be very limited and we will subsequently show that this is indeed the case. 



4 Model with Dynamic Risk Aversion 

4.1 Heterogeneous and dynamic risk aversion 

To extend the baseline model, we allow agents to have heterogeneous risk aversion indices (de- 
grees), which vary with time. This reflects the fact that in a real financial market investors 
have different risk attitudes and the investors' sentiment change with time. We assume that the 
risk aversion index of each agent follows an independent bounded random walk with a constant 
variance <5 2 : 

7M = Ji,t-l + 8z itt , Ji,t G [70, 7u] ( 25 ) 

where z^t is an i.i.d. Gaussian variable with mean zero and unit variance for agent i, 7o(> 0) is 
the lower boundary, and 7«(> 70) the upper boundary. It's easy to relate the value of the index 
at time t + r to the value at time t: 

T 
t=l 

where Si >T = Ylt=i z i,t ^ s the change of risk aversion index of agent i from time t to time t+r, which 
can be either positive or negative. It's worthwhile to note that in real markets, the dynamics of 
investors' risk aversion attitudes may be more complicated than a simple random walk process 
we assume here. However, simplifying and idealizing of the real situation helps us to stay focused 
on the main purpose of investigating the impact of investors' fluctuating risk aversion on the 
price dynamics. 

4.2 Price setting equation with dynamic risk aversion 

Upon substitution of Eqn. (|26[) into Eqn. Q, we have: 

v^N 1 ( E itt + T (pt+T+i+dt+T+i) X \ 

_ ^ (7M+agi.r)l ^ t+T (l+r) ^7 (07 , 

Pt+T ~ 1 / i?i, t (pt+i+rft+i) 1) Pt { > 

Comparing Ean.(|27|) to Eqn. (|2*T|) (7^4 = 7i i t+ r = 7«) for the case of fixed constant risk aversion, 
we see that there is an extra term, SSi tT , in the price setting equation in the case of DRA. Since 
Sj jT can be either positive or negative and its value changes with time, 7^ + SS^ T deviates from 
7i t and fluctuates with time. This fluctuating weighting factor (representing agent's fluctuating 
risk aversion) acts like an 'amplifier" of the price deviation induced by the error in agents' price 
estimation, and therefore results in excess price fluctuation. \S^ T \ ~ iV(0, y/r), for yfr5 S> 1, 
li,t+r{li,t + 8Si jT ) and 7^ can differ substantially, resulting in a large deviation of pt+ T from 
Pt- Our numerical results, presented in the next session, clearly show that it's this risk aversion 
dynamics that gives rise to the excessive price fluctuations and the associated volatility clustering. 
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4.3 The range of DRA indices 

We now examine the range of possible relative risk aversion indices. The choice of the range is 
important for modeling investors' decision-making; it has big impact on the price dynamics, as it 
directly affects investors' demand of the risky asset. The lower the index, the less risk-averse the 
investor is (thus the higher the demand of risky asset); and vice versa. Thus the risk-aversion 
attitude has great impact on the price dynamics through its influence on the demand. Some 
empirical and experimental studies reported that for a "typical" investor, the value of the risk- 
aversion index 7 is in the range of 0-2 (Mehra and Prescott (1985), Friend and Blume (1975), Levy, 
Levy and Solomon (2000)). Mehra and Prescott (1985) used a value of risk aversion index with 
an upper limit of 10 in their treatment of the issue of the "Equity Premium Puzzle" . However, 
to "explain" the "Equity Premium Puzzle" of NYSE over 50 years of U.S. postwar period, one 
needs a relative risk aversion index of 250 if a consumption-based model is used(Cochrane (2005))! 
These empirical results show that it is better to model the risk aversion with a range of indices, 
instead of a fixed value. The range we specified consists of an upper bound and a lower bound 
for the random walk describing DRA indices. 

5 The Simulation Results and Analysis 

5.1 The setup 

In our simulation we choose the number of agents N = 100, the number of predictors each agent 
has, M = 2. Setting different number of agents produces similar results. The initial risk aversion 
indices 7^0 are all set to 1.0 for the baseline model and are set to 7^0 & [0.2, 4] for the model with 
DRA. The bounds for the index of DRA are 7^ £ [10 -5 , 20], the risk-free interest rate is r = 5%, 
the dividend growth rate is = 2%, and the weighting coefficients for the variance of estimation 
is 6 = 1/250. The lags used in the price estimators are Ljj G [2,250], and we set <j p+ d=l% for 
all agents. 

5.2 Simulation price and trading volume 

Let's first take a look at how excess price fluctuations emerge from a dynamic risk aversion 
process. Fig.l shows the simulation time series of the price and the trading volume generated 
from the model with fixed constant risk aversion (CRA) and the model with dynamic risk aversion 
(DRA). 

From the figure we see clearly that the DRA leads to increased fluctuations in both the price 
and the trading volume. To have both qualitative and quantitative picture of the impact of the 
DRA on the price dynamics, we examine the key stylized facts in the following subsections. 

5.3 Autocorrelation function 

One of the stylized facts observed in real financial data is that their autocorrelation functions 
(ACF) usually start with a low value (from pi) and decay very slowly with increase of time step 
for the squared or absolute-valued returns. For an almost-Gaussian process, the values of its 
ACFs for absolute- valued return are close to zero and independent of the time steps. In Fig. 2, 
we compare the ACFs for absolute-valued returns for the series generated by our baseline model 
(with CRA) and the DRA model, the series of real data (DJIA and SP500 Index), and Gaussian 
process. 
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The figure clearly shows that while the ACFs generated from our baseline model (CRA) is 
close to that of Gaussian process (Gauss), the results generated from our DRA model are very 
close to the real financial data (DJIA, SP500). 

5.4 Excess volatility 

The second key stylized facts we examine is the excess volatility (or fat-tails) of returns, which 
measures the price fluctuation of real financial series. In Fig.3 we plot the distributions of returns 
(in different time-steps) from our baseline model and the DRA model with the parameters set 
according to a normal market condition. For comparison, we also plot the return distribution of 
DJIA and the return distribution generated by a simple Gaussian process. 

These plots show that, for all different time periods, the return distributions from our baseline 
model are very close to those of the Gaussian process; in contrast, the results from the model with 
DRA are close to the real DJIA data. In the context of our model, it is clear that the dynamic 
risk aversion leads to excess volatility or a fat-tail in the return distribution, which is one of the 
most important characterization of real financial time series (Mandelbrot (1963), Fama (1965)). 

To further examine the fat tail of the distribution, we plot, in Fig. 4, the Kurtosis as a function 
of the squre root of variance of risk aversion, 5. From these plots, we see that the risk aversion 
dynamics can change the return distribution significantly from a Gaussian distribution (which has 
K=3.0). In addition, the smaller the lag r, the larger the Kurtosis generated; this is consistent 
with the empirical observations in real financial data. These values of the Kurtosis, together with 
the standard deviation and skewness are listed in Table 1. Note that the statistics generated from 
our DRA model are quite close to that from DJIA daily data. 

5.5 Volatility Clustering 

Volatility clustering is another important characteristics of financial time series. Here we use 
the conditional probability measure developed recently by Chen et al. (2005) to examine the 
volatility clustering. The method uses the return distribution conditional on the absolute return 
in the previous period to describe a functional relation between the variance of the current return 
and the absolute return in the previous period. If the volatility in asset returns is clustered, it 
will be proportional to the volatility in the previous period, the proportionality constant reflects 
the strength of volatility clustering. 

Fig. 5 shows how the current volatility depends on the volatility of the previous period for a 
random walk process, the baseline market model, DRA model, and DJIA daily data. From the 
figures we see clearly that the baseline market model with fixed constant risk aversion (CRA) 
produces very low volatility clustering (the curve is flat if there is no volatility clustering, such 
as the case with the random walk model (Gauss)). In contrast, the volatility clustering from the 
DRA model is significantly higher, and it is very close to the one generated from DJIA daily 
prices. The plots also show that, for our DRA model, the smaller the time period (step) of the 
return, the stronger the volatility clustering; and vice versa; this is consistent with the empirical 
observations of real market data. 

6 The SFI market model with dynamic risk aversion 
6.1 Brief introduction to SFI market 

To test the impact of DRA on other baseline model, we use the well-known Santa Fe model of 
artificial market (Arthur et al. (1997), LeBaron, Arthur and Palmer (1999)). We first give a 
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brief summary of the SFI market below. In the SFI market a constant absolute risk aversion 
(CARA) utility function (U(ct, 7) = — e c * 7 ) is assumed for each agent, the demand and price 
setting equations in fact have similar forms as in Eqn. (J2j) and ®. 
The dividend process is assumed to be an AR(1) process: 

d t = d + p{d t -d) + s t (28) 

where ej is an i.i.d. Gaussian with zero mean and variance a e 

The SFI market assumes that each of the N (=25) agents at any time possesses M (=100) 
linear predictors and uses those that best matches the current market state and have recently 
proved most accurate. Each predictor is a linear regressor of the previous price and dividend, 
E{p t j r \ + dt+\) = a(pt + dt) + b\ it uses a market state "recognizer" vector consisting of J (=12) 
elements, each taking a value of either 0, 1 or #(match any market states). The market status 
is described by a state vector consisting of J binary elements, each taking value of either 1 (its 
specified market condition exists) or (the market condition does not exist). The elements of 
market state can represent any important market discriminative information, including macro- 
/micro- economic environment, summary of fundamentals, and market temporal trends, etc. At 
each time step, only those predictors which match all their J elements to the corresponding J 
elements of market status are eligible to be used and are called "active" predictors. 

The variance of estimation a\ t for each agent is assumed to update with a moving average of 
squared forecast error, defined in Eqn. (|17j). 

Agents learn to improve their performance by discarding the worst (20%) predictors and 
developing new predictors via a genetic algorithm. This ensures the market some dynamics. 

In the SFM, the market conditions are specified as: 

1-6 elements represent "current price x interest rate / dividend > 0.25, 0.5, 0.75, 0.875, 1.0, 
1.125. 

7-10 elements describe "Current price > 5-period moving average of past prices (5-period 
MA), 10-period MA, 100-period MA, 500-period MA. 
11th element always 1; 
12th element always 0; 

The regressor's parameters {a, b} are set to be randomly and uniformly distributed within 
the ranges: a G (0.7, 1.2) and b G (-10, 19.002). 

The risk-free interest rate is set to 5%, and for the dividend process: the auto-regression 
coefficient p = 0.95, d = 10, and a e = 0.0745. 

The weighting coefficients for the variance of estimation, = 1/75, 1/150 for faster and slower 
learning respectively. The genetic algorithm is invoked (on average) every T e = 250 periods 
(faster learning) and 1000 periods (slower learning). For more detailed justification for choosing 
the parameter values, see LeBaron et al. (1999). 

With these setups, we have checked that the simulation stock price time series and its statis- 
tical properties generated are in fact similar to those of our baseline model. 

6.2 Numerical results of SFI market model with DRA 

The dynamics of risk aversion can be similarly incorporated into the SFI market model. Fig. 6 
plots the simulation price time series for different DRA variance 5 2 . These plots show clearly the 
impact on price fluctuations from the DRA. To have a quantitative picture on how the excess 
volatility emerges from the DRA, we plot in Fig. 7 the functional relation of the Kurtosis vs 
variable 5. The results are very similar to those plotted in Fig. 4, which were obtained with our 
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much simpler baseline model. 

We have checked other key features and found that the SFI-DRA model gives the similar 
results as our DRA model which is based on a much simpler baseline model. This suggests that 
DRA is the key mechanism for the emergence of the key stylized facts, and the impact of DRA 
does not depend on the structures of the baseline models. Therefore the price impact of investors' 
DRA we have studied is generic. 

7 Summary 

We have presented a simple multi-agent model of a financial market which incorporates the dy- 
namics of risk aversions. We assume that the index of DRA follows a simple independent bounded 
random walk with a constant variance 5 2 . We demonstrate that such dynamics is directly re- 
sponsible for excess volatility and the associated volatility clustering. We compare the numerical 
results from our model with the results obtained by analyzing the DJIA daily data and show that 
the simulation data reproduce most of the "stylized" facts, such as excess volatility (measured by 
fat tail and high peak of return distribution) , volatility clustering measured by conditional return 
distribution. We have also tested the DRA on the Santa Fe market model and obtain similar 
results. This suggests that the impact of DRA among heterogeneous agents we introduced here 
does not depend on the structure of the particular baseline model used. The degree of excess 
volatility is essentially controlled by the parameter 5. Thus S can be used as a key market sen- 
timent parameter, in conjunction with the other market indicators such as average return r and 
the average volatility o~o, to characterize the financial market. We hope that our results presented 
here will provide new insights into the dynamics of asset price fluctuations governed by investors' 
fluctuating sentiments. 
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Table 1: Some statistics for DRA model ( 5=0. 01) and DJIA 



Lag (r) 


S.D. 


Skewness 


Kurtosis 


Model S = 0.01 


l 


0.015004 


-0.454483 


48.406593 


2 


0.018313 


-0.372956 


35.233093 


4 


0.021685 


-0.460952 


26.051126 


8 


0.026179 


-0.492277 


17.309658 


16 


0.031840 


-0.424137 


11.376849 


32 


0.039121 


-0.327192 


7.444087 


DJIA 


l 


0.010876 


-1.190497 


40.2330 


2 


0.015710 


-1.100639 


29.9739 


4 


0.022213 


-0.950496 


18.9213 


8 


0.032103 


-0.989078 


13.6606 


16 


0.046404 


-0.995792 


12.0792 


32 


0.066602 


-0.706621 


9.5892 
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Simulation price time series for different 5 (from bottom) = O, 0.005, 0.01, 0.015, and 0.02 




Simulation trading volume time series for different 5 =0 (bottom) and = 0.1 (top) 




Figure 1: The time series of price and trading volume from the models with constant (S = 0) and dynamic 
(5 ^ 0) risk aversion. For the sake of clarity, the time series with different <5s were vertically shifted, e.g., the 
trading volume for 8 =0.01 was upward shifted by 50. 




200 250 300 

time step t 



Figure 2: ACFs of different time steps for return series of: the baseline model (CRA), the DRA model, DJIA 
Index, SP500 Index, and random walk process (Gauss). 
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Figure 3: Return distributions for the model with constant risk aversion (CRA) (5=0) and the DRA model 
(5=0.01), Gaussian process (Gauss) and real data of DJIA. 
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Figure 4: Kurtosis of simulation price series for different variables 8 of the DRA index. 
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Figure 5: The volatility clustering measured by the standard deviation of current return vs the absolute return 
of previous period. 




Figure 6: Simulation price series for different variances of DRA processes with Santa Fe model. The model 
parameters are: 9 = 1/75, T e = 250. The prices for different 8s have been vertically shifted to make the comparison 
clearer. From the bottom to top, the price series are respectively for 8 — 0, 0.01, 0.02, 0.03, 0.04 and 0.05. 
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